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Abstract 

We study the behaviors of entanglement entropy and vacuum expectation 
value of Wilson loop in the Stiickelberg holographic insulator/superconductor 
model. This model has rich phase structures depending on model parameters. 
Both the entanglement entropy for a strip geometry and the heavy quark poten- 
tial from the Wilson loop show that there exists a "confinement / deconfinement" 
phase transition. In addition, we find that the non-monotonic behavior of the 
entanglement entropy with respect to chemical potential is universal in this 
model. The pseudo potential from the spatial Wilson loop also has a simi- 
lar non-monotonic behavior. It turns out that the entanglement entropy and 
Wilson loop are good probes to study the properties of the holographic super- 
conductor phase transition. 
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1 Introduction 

The entanglement entropy is expected to be a key quantity to understand some properties 
in quantum field theories and in many-body physics (see, for example, Refs. [1, 2]). For a 
given system, the entanglement entropy of one subsystem with its complement is defined 
as the von Neumann entropy. The entanglement entropy of the subsystem measures how 
the subsystem and its complement are correlated each other. As a von Neumann entropy, 
the entanglement entropy is also directly related to the degrees of freedom of the system. 
In addition, in quantum many-body physics, the entanglement entropy is a good quantity 
to characterize different phases and associated phase transitions. However, the calculation 
of entanglement entropy is very difficult except for the case in 1 + 1 dimensions. In the 
spirit of the AdS/CFT correspondence [3, 4, 5, 6], a holographic proposal to calculate the 
entanglement entropy for strongly coupled field theories has been presented in Ref. [7] (for 
reviews see [8, 9]). 

On the other hand, while various aspects of holographic superconductor models have 
been intensively studied (see, for example, Refs. [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 
20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31]), the study of entanglement entropy in the 
holographic superconducting phase transition is just initialed. Refs. [32, 33] studied the 
behaviors of entanglement entropy for a strip geometry in the holographic s-wave and p- 
wave conductor /superconductor models. It turns out that the entanglement entropy is 
a good probe to investigate the holographic phase transition, the entanglement entropy 
behaves like the thermal entropy of background black holes, it can indicate not only the 
appearance, but also the order of the phase transition. In a recent paper [34], the authors 
investigated the behavior of the entanglement entropy for a strip geometry in a simple 
holographic insulator /superconductor model at zero temperature. In this model, the phase 
transition is a second order one. It was found that the entanglement entropy as a function 
of chemical potential is not monotonic in the superconducting phase. The entanglement 
entropy at first increases and reaches its maximum at a certain chemical potential and then 
decreases monotonically as chemical potential increases. The non-monotonic behavior of 
entanglement entropy in the superconducting phase looks strange. Due to the lack of 
details of dual field theory, we did not give a convinced interpretation, although some 
possible causes on this behavior were discussed there [34]. 
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In the holographic insulator/superconductor model [35], the normal insulator phase is 
described by a pure AdS soliton solution. As one increases the chemical potential, the 
pure AdS soliton background will become unstable to develop some kind of scalar "hair" . 
The emergence of scalar "hair" induces the symmetry breaking and gives a finite vacuum 
expectation value of the dual operator in the field theory side. The soliton solution with 
scalar "hair" describes a superconducting phase. The scalar "hair" therefore plays the role 
of the order parameter in the holographic phase transition. 

In this paper we are going to further study the behavior of the entanglement entropy in 
the holographic insulator /superconductor phase transition by generalizing the discussion 
to the case in the Stiickelberg holographic insulator /superconductor model [36]. The aim 
is two-fold. On the one hand, we would like to see whether the non-monotonic behavior is 
universal or not. On the other hand, the phase structure is rather rich in the Stiickelberg 
holographic insulator/superconductor model. As we will see shortly, depending on the 
model parameters, the superconducting phase transition could be second order or first 
order, and further a first order phase transition also occurs in the superconducting phase. 

The Stiickelberg holographic insulator/superconductor model is labeled by two param- 
eters. One is /3 determining the strength of the back reaction. The other is ( determining 
the form of interaction between the scalar field and Maxwell field. Both of them can change 
the order of the phase transition. As £ is vanishing, an equivalent model has been studied 
in Ref. [37]. It shows that as one increases the strength of the back reaction, the order 
of the phase transition is changed from second order to first order. For the intermediate 
strength of the back reaction, although the insulator /superconductor transition is second 
order, a new phase transition emerges in superconducting phase. We fix the strength of 
the back reaction and study the effect with the change of (. Interestingly, the resulting 
phenomenon is qualitatively similar to the above case. The order of the phase transition 
is second order for small £ and first order for large (. For the intermediate value of £, the 
phase transition is second order, however, the grand potential in superconducting phase 
develops a "swallow tail" indicating a new phase transition. 

We calculate the entanglement entropy for a strip geometry in this model, and find that 
the entanglement entropy can indicate not only the appearance of phase transition, but also 
the order of the phase transition. Further, no matter the order of the phase transition, the 
entanglement entropy versus chemical potential is always non-monotonic in the supercon- 
ducting phase. More precisely, at the beginning of the transition, the entropy increases and 
reaches its maximum at a certain chemical potential and then decreases monotonically. It 
indicates that the non-monotonic behavior of entanglement entropy as a function of chem- 
ical potential is universal in the holographic s-wave insulator/superconductor model. 

Apart from the entanglement entropy, there is another nonlocal quantity, Wilson loop, 
which can describe some properties of gauge field theories. In the AdS/CFT correspon- 
dence, the vacuum expectation value (VEV) of Wilson loop can also be calculated holo- 
graphically. In order to give further insights into the holographic insulator/superconudctor 
phase transition, we study the behaviors of temporal Wilson loop and spatial Wilson loop 
across the holographic phase transition. From the VEV of Wilson loops, we extract the 
heavy quark potential which describes the interaction strength between quark and anti- 
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quark. From the entanglement entropy and the heavy quark potential, it shows that there 
exists a "confinement /deconfmement" phase transition in the dual field theory describing 
the insulator /superconductor model. In addition, different from the phenomenon observed 
in the entanglement entropy, the (pseudo) heavy quark potential versus chemical potential 
in the superconducting phase will show the monotonic behavior or non-monotonic behavior, 
depending on the model parameters. 

This paper is organized as follows. In section (2), we introduce the Stiickelberg holo- 
graphic model. In section (3), the insulator/superconductor phase transition and its ther- 
modynamics are investigated in detail. Section (4) is devoted to investigating the behavior 
of entanglement entropy in the holographic model. In section(5), we study the behavior 
of temporal Wilson loop and spatial Wilson loop with respect to chemical potential and 
distance between quark and antiquark. The conclusion and some discussions are included 
in section (6). 



2 Stiickelberg Insulator/Superconductor Model 

The Stiickelberg holographic superconductor model which contains a real scalar field ip, a 
real pseudoscalar field p and a Maxwell gauge field reads [36, 38, 39] 

/I 12 1 

d 5 Xy/^g[— (K + — ) + -^Lmatter], 
Zk L g /j\ 

where k is the gravitational constant and J 7 is a function of ip. In this paper we consider 
a simple case J r (^) = ip 2 + (ip 6 , where ( is a model parameter determining the interaction 
form and is assumed to be non-negative to ensure the positivity of the kinetic term for p. 
For other forms of J 7 , see [36, 38, 39]. The local U(l) gauge symmetry in this theory is 
given by 

p^p + 9(x»), A^A^ + V^). (2) 

We use this symmetry to set p = from now on. Here we define a parameter = n/g 
which measures the strength of the back reaction of the matter fields on the background 
geometry. 

The equations of motion coming from the above action are 

- —g^ - —g^y = (3 2 [F^F V X + 2V M ^V^ + 2(^ 2 + C^)^A V + g^L matter ], 

V U F U » = 2(</> 2 + C^ 6 )^. 
Our ansatz for the metric and matter fields are given by 

T2 J 2 

ds 2 = —— + r 2 (-f(r)dt 2 + dx 2 + dy 2 + gMe^dr) 2 ), (4) 
r 2 g(r) 
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vj) = i){r), A = 4>(r)dt, (5) 

where g(r) vanishes at the tip of the soliton r = r . Further, in order to avoid a conical 
singularity at the tip tq, rj should be made with an identification 

X(rp) 



ri~ v + r, r = 2 —_ - s . (6) 



4nLe~ 

The independent equations of motion in terms of the above ansatz are deduced as follows. 
/// ,/ 5 , /' , 9' X\,, , L 2 4> 2 - L 2 m 2 

0" + (--^ + ^-^)0'-^(^ 2 + ^ 6 ) = o, 

r 2j g 2 r z g 

/" + (-- £ + + + 4/3V 2 )/ - 2 -^f = 0, 
r 2j 2 r r z 

,3 /' ,12 /Y oo2(/2 /?V\ 2L 2 /3 2 2 2 ^ 2L 2 /3 2 mV-12 n 

r 2] r z 2j r z j r^j r z 

,3f , £l . , ,3 , , 4L 2 /3 



2/32 J2 



+ /v - (- + -)f + ^f^c^ + c^ 6 ) + w 2 fr 



r r g ' ' r A g 



(7) 



In our numerical calculations, we choose m 2 L 2 = — ^ and work in units with L = 1. How- 
ever, this analysis can be generalized to any mass case above the Breitenloher-Freedman 
bound m 2 L 2 > —4. In order to match the asymptotical AdS boundary, the matter and 
metric fields near the boundary r — > oo should have the form 

V = + ~TH + •••' ^ = ^0 n + • • • , 

/ 1 i ■/ 4 i 1 i 5,4 i ^ 4 1 

/ = i + ^ + ..., <7 = i+^ + ..., x = ^ + ..., 

where ipo, <po, <p2, fi, #4 and \i are all constants. It is well-known that in five- 
dimensional AdS space-time, when —4 < m 2 L 2 < —3, the scalar field admits two dif- 
ferent quantizations related by a Legendre transformation. ip can either be identified as 
a source or an expectation value. In this paper we consider it as source and set ipo = to 
accomplish spontaneous symmetry breaking. The quantities <f)Q, <p2 and ipi are related to 
the chemical potential /1, charge density p and the vacuum expectation value of the scalar 
operator O which has scaling dimension A = | in the dual field theory on the boundary, 

i.e., „ = 0o, P = g0 2 and (O) = = 



3 Thermodynamics and Phase Transition 

In gauge/gravity duality the grand potential (Gibbs free energy) Q of the boundary thermal 
state is identified with temperature times the on-shell bulk Euclidean action. Namely 
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f2 = T S Euclidean- The Euclidean action must include the Gibbons-Hawking boundary term 
for a well-defined Dirichlet variational principle and further a surface counter term for 
removing divergence 

/I 12 1 If 
d5x V9[7nO Z + 71 ) + -^ L matter] J / d^xVhK + S ct , (9) 

where h is the induced metric on the boundary r — > oo, and /C is the trace of the extrinsic 
curvature. 5 c t is the counter term given by 

Set = - 2 [ d*xVh(j + /? 2 i^ 2 ). (10) 

For our soliton solution (4), there is no horizon and associated Hawking temperature and 
thermal entropy vanish. But for the Euclidean sector of the solution (4), one can introduce 
an arbitrary inverse temperature (1/T) as the period of the Euclidean time coordinate. 
Due to the fact that the soliton solution is static, the integration over the Euclidean time 
in the Euclidean action gives an inverse temperature factor 1/T, which just cancels the 
temperature factor in the grand potential and leads a finite grand potential. Using the 
on-shell condition and the expansion of matter and metric functions at infinity r — > oo, 
the grand potential Q is found to be given as 

rr = 94, (11) 

V3 

with V3 = J dxdydr]. Since we have scaled T to be irL, one has g± = —1 in the normal 
insulator phase. 

Figure (1) shows the behavior of the condensate with respect to chemical potential for 
the back reaction parameter /3 = 0.25 1 . When ( is small, we see from Figure (1) that as the 
chemical potential \i exceeds a critical value \i c the condensate emerges, which implies a 
superconducting phase appears. On the other hand, when \i < fi c , the scalar hair vanishes. 
This is identified as the insulator phase, since the system has a mass gap ~ 1/T. The 
insulator /superconductor phase transition here is typically second order and in this case, 
the critical chemical potential /x c does not depend on (, the critical chemical potential is 
given by /i c ~ 1.888 2 . 

As we increase ( to £ ~ 0.5, an additional complication appears. The condensate versus 
chemical potential does not have a monotonic behavior (see the red curve in Figure (1)). 
As we can see clearly in the left plot of Figure (2), the grand potential with respect to 
chemical potential develops a "swallow tail", which is a typical signal for a first order 
phase transition. It implies that there is a new phase transition within the superconduct- 
ing phase. If we continue to increase £, the insulator/superconductor phase transition will 



1 In our numerical calculations of this paper we always take /3 = 0.25 and show rich phase structure by 
changing the parameter £. 

2 When C is beyond a certain value, the superconducting phase transition will become a first order one. 
In that case, the critical chemical will depend on £. 
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2.0 2.5 3.0 3.5 

Figure 1: The condensate of operator O versus chemical potential \x for /3 = 0.25. Lines 
from right to left are for £ = (purple), £ = 0.2 (green), £ = 0.65 (red) and ( = 1.6 (blue), 
respectively. T is scaled to be tiL. 

become first order when £ > 1, which can be seen in the right plot of Figure (2). 3 This is 
qualitatively similar to the case in Ref. [37] where the effect of back reaction on the insu- 
lator/superconductor was studied in the case with vanishing (. The condensate becomes 
larger with the increase of £, which means that the scalar hair will be more difficult to be 
formed in the AdS soliton background with larger (. 



4 Entanglement Entropy 

Now let us begin to study the behavior of entanglement entropy in this holographic super- 
conductor model. The holographic method to calculate entanglement entropy is as follows. 
Consider a strongly coupled field theory with gravity dual, the entanglement entropy of 
subsystem A with its complement is given by searching for the minimal area surface 7^ 
extended into the bulk with the same boundary dA of A. Then the entanglement entropy 
of A with its complement is given by the "area law" [7] 

S A =^Area( lA ), (12) 

where k is the gravitational constant. 

We consider the subsystem A with a straight strip geometry with a finite width i along 
the x direction, along the rj direction with a period T, but infinitely extending along the y 

3 It is worth pointing out here that if the back reaction parameter f3 is larger enough, the transition will 
always be first order, no matter how small the value of £ is chosen. However, we do not consider the case 
in this paper. 
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Figure 2: The grand potentials of the soliton with scalar hair (solid) and the soliton without 
scalar hair (dashed black) with respect to chemical potential for £ = 0.65 (left plot) and 
( = 1.6 (right plot), respectively. T is scaled to be nL and = 0.25. The left plot has a 
second order phase transition at /x ~ 1.888, but a discontinuity at fi ~ 1.962 (denoted by 
dotted vertical green line) in superconducting phase. The right plot shows the typical first 
order phase transition and the critical chemical potential denoted by vertical dashed pink 
line is fi c ~ 1.813. 



direction. The subsystem A sites on the slice r = -, where e — > is the UV cutoff. The 
holographic dual surface 7^ is defined as a three-dimensional surface 

* = 0, r = r{x), ~<y<5 (R-Kx,), < 77 < T. (13) 

We are first interested in the case that the surface is smooth. The holographic surface 7^ 
starts from x — | at r = -, extends into the bulk until it reaches r = r*, then returns back 
to the AdS boundary r = -ata; = — |. In this case, the entanglement entropy is given by 

ST n ect = ^L RT f r^e-^ ^ = 2nL R 1 ^ co 

A « 2 Jr. v /r 6 ^(r)e-xM _ r 6 ,g{r,)e~^ k 2 e 2 

where the UV cutoff 1/e has been taken into consideration. The width t of the subsystem 
A and r* are related by 

t\ T 

=dr. (15) 



r 2 / n ( r \( r6g(r)e-x(r) 

In fact, there are two solutions for the connected configuration. See Figure (3). In addition, 
there is also a disconnected configuration describing two separated surfaces located at 
x = ±| respectively and extending to the bulk and reaching at the tip of the soliton 
geometry. The entropy for this disconnected geometry is independent of £, and given by 

^connect = ^L RV [' re ~^ dr = ^ RT (1- + S dlSCOn ) . (16) 
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Figure 3: The entanglement entropy as a function of strip width £ for /3 = 0.25, £ = 0.65 and 
\i ~ 2.276. The dashed purple and solid blue curves come from the connected configuration, 
while the solid red one comes from the disconnected configuration. The lowest curve is 
physically favored compared with others. In this figure £ cr u — 0.66L and £ m ax — 0.71L. 

We find that the entanglement entropy with respect to strip width I behaves quite 
similar for different choice of parameters, i.e., f3, ( and /i. We draw Figure (3) with f3 = 0.25, 
C = 0.65 and fi ~ 2.276 as a concrete example. The connected configuration does not exist 
when i > imax- The physical solution comes from the trivial disconnected geometry and the 
entropy is independent of strip width. On the contrary, there are three different branches 
when £ < £ max , i-e., the upper branch (dashed purple), the lower branch (solid blue) and the 
middle branch (solid red). The first two curves correspond to the connected surface, while 
the third curve comes from the disconnected one. We trace over the physical entropy by 
always choosing the lowest branch. As we can see from Figure (3), the lower branch (solid 
blue) is finally favored as we decrease I. Hence, there is a critical value £ cr u below which 
the lower branch is physically favored. Thus, as we change £, a phase transition occurs at 
£ cr it, which is just the so called "confinement/deconfinement" phase transition [40, 41, 42]. 
To be precise, for £ < £ cr u, the entanglement entropy comes from the connected surface 
and exhibits non-trivial dependence on £, which describes a "deconfiement" phase. For 
£ > £crit, the entropy is dominated by the disconnected configuration and is £ independent, 
which indicates a "confinement" phase. 

To summarize, there are totally four "phases" in the dual boundary field theory, i.e., the 
insulator phase, superconductor phase, and their corresponding confinement/deconfinement 
phases. These phases are characterized by the chemical potential fi and strip width £. In 
particular, the strip width controls the "confinement/deconfinement" phase transition 4 . 

4 Strictly speaking, the term phase transition here is inappropriate since the system itself, i.e., the state 
of the boundary field theory, does not change at all as one changes i. However, this observed behavior 
here is quite similar to the one in the thermodynamic confinement/deconfinement phase transitions and 
therefore we follow Rcfs. [40, 41, 42] and adopt the terminology "phase transition" to describe this behavior. 




max 
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Figure 4: The entanglement entropy as a function of chemical potential for £ = (purple 
curves) and ( = 0.2 (green curves) at fixed £/L = 0.5 (within "deconfining phase", left 
plot) and £/L — > oo (within "confining phase", right plot), respectively. The transition 
from insulator to superconductor is second order here. In both two plots, the solid curves 
come from superconducting phase, while the dashed black lines come from insulator phase. 

It is instructive to study how the entanglement entropy changes with chemical potential 
by fixing strip width i. We first focus on the second order phase transition case, which 
is presented in Figure (4). We can see that, at the beginning of the phase transition, the 
entanglement entropy increases continuously with chemical potential both in "deconfine- 
ment phase" and "confinement phase" and reaches its maximum at a certain chemical 
potential fi m ax, then it decrease monotonically. Furthermore, the entanglement entropy is 
continuous at critical chemical potential fi c , but its slop has a discontinuous change at \x c . 
The behavior of the entanglement entropy across the phase transition point indicates that 
the transition is a second order one. 

When 0.5 < ( < 1, a jump in the condensate (the red curve in Figure (1)) and a 
"swallow tail" in the grand potential (see the left plot in Figure (2)) appear, although the 
phase transition is still a second order one at /x c . Interestingly, we can also see a jump 
in Figure (5) which shows the entanglement entropy as a function of chemical potential 
at fixed strip width. Tracing the physical curve, we find that the entropy increases at 
the beginning of the transition, then there is a sudden jump at the chemical potential 
/io — 1.962. The sudden jump of entropy indicates a first order phase transition there. 

For sufficiently large (, the insulator /superconductor transition becomes first order. 
The entanglement entropy with respect to chemical potential is presented in Figure (6). 
Comparing with Figure (4) and (5), we can see a dramatic change in the first order case. 
Although the entropy in superconducting phase also behaves non-monotonically with re- 
spect to /i, the entropy as well as its slop at the critical point have a discontinuous jump. 
Once again, it is a signal of a first order transition. 

In Figure (7), we plot the critical length £ crit of the "confinement /deconfinement" phase 
transition with respect to the chemical potential in the superconducting phase. For small 
(, £ cr it first increases and forms a peak, then it will increase continuously for large chemical 
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Figure 5: The entanglement entropy as a function of chemical potential for ( = 0.65 at fixed 
i/L = 0.5 (within "deconfinement phase", left plot) and i/L — > oo (within "confinement 
phase", right plot), respectively. The transition from insulator to superconductor is second 
order, but a discontinuity at fiQ ~ 1.962 (denoted by dotted vertical line) within the 
superconducting phase. In both two plots, the solid red curves are from superconducting 
phase, while the dashed black lines are from insulator phase. Trace the physical curve by 
choosing the dashed black line for [i < fi c , then choosing the red curve which has a jump 
from point A to point B at /i — 1.962. 

potential. As we increase (, the position of the peak moves toward to the region with small 
chemical potential and will be finally cut off since the transition is first order for large £, as 
we can see that the dotted part of the blue curve in Figure (7) is not thermodynamically 
favored. Here an remarkable point is that the critical length as a function of the chemical 
potential is not monotonia This is of course related to the fact that the entanglement 
entropy is not a monotonic function of the chemical potential. 

We can see from Figures (4), (5) and (6) that although the phase structures are different 
in three cases, the entanglement entropy with respect to chemical potential are always non- 
monotonic in superconducting phase. The discontinuity or jump at critical point indicates 
some kind of significant reorganization of the degrees of freedom of the system, since new 
degrees of freedom are expected to emerge in new phase. The non-monotonic behavior 
of the entanglement entropy is interesting. To further study its implication, in the next 
section, we will discuss the behavior of another non-local quantity, Wilson loop, in the 
superconducting phase. 

5 Wilson Loop 

From the VEV of Wilson loop one can extract the interaction potential between a quark- 
antiquark pair in gauge field theory. Therefore the VEV of Wilson loop is thought of as 
a proper quantity to describe the confinement /deconfinement phase transition of gauge 
field theory. In the AdS/CFT correspondence, a proposal for the VEV of Wilson loop is 
presented in Refs. [43, 44, 45]. In this section, following the proposal we will discuss the 
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Figure 6: The entanglement entropy as a function of chemical potential for ( = 1.6 at fixed 
£/L = 0.5 (within "deconfinement phase", left plot) and £/L — > oo (within "confinement 
phase", right plot), respectively. The transition is first order and the critical chemical 
potential denoted by vertical dashed pink line is fi c ~ 1.813. In both two plots, the 
solid blue curves are from superconducting phase, while the dashed black lines are from 
insulator phase. Trace the physical curve by choosing the dashed black line for fi < fi c , 
and for fi > fi c , choosing the blue curve starting from point D. 

behavior of the VEV of Wilson loop in the holographic insulator /superconductor phase 
transition. We will consider two kinds of Wilson loop: One is a temporal Wilson loop; the 
other is a spatial Wilson loop. In both cases, we will show the behavior of the interaction 
potential with respect to the distance £ between quark and antiquark and to the background 
chemical potential /x. 

5.1 Temporal Wilson Loop 

In SU(N) gauge theory, the interaction potential for heavy quark-antiquark (QQ) pair can 
be calculated from the Wilson loop 

W[C) = ^-TrPexp[i I A^daf], (17) 
N Jc 

where is the gauge field, the trace is over the fundamental representation, P stands 
for path ordering. C denotes a closed loop in space-time, which is a rectangle with one 
direction along the temporal direction of length T and the other spatial direction of length 
I. The Wilson loop describes the creation of a QQ pair with distance £ at some time to = 
and the annihilation of this pair at time t = T. For T — > oo, the VEV of the Wilson loop 
goes as (W(C)) oc e~ TV QQ. In terms of the AdS/CFT dictionary, the VEV of the Wilson 
loop in four dimensions should be equal to the string partition function on the curved 
space, with the string world sheet ending on the contour C at the boundary of the curved 
space [43] 

(W 4d [C\) = Z™[C] ~ , (18) 
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Figure 7: The critical length as a function of chemical potential in superconducting phase, 
where the purple curve (up) is for ( = 0, the green curve (middle) is for ( = 0.2 and 
the blue curve (down) is for ( = 1.6. The dotted part of the line for ( = 1.6 (blue) is 
not thermodynamically favored, trace the physical curve by always choosing the solid line. 
Note that the curve for ( = (purple) will continue to increases large chemical potential 
as other curves. 



where S x -t is the classical world sheet action of the Nambu-Goto form 

S x - t = J drdcr^DetXab, (19) 

and a' is the string tension with dimension [energy]~ 2 , and Xab is the induced world sheet 
metric with a, b the indices in the (r = t, a — x) coordinates on the world sheet. In this 
subsection, we follow the standard procedure [43, 44, 45] to extract the static heavy quark 
potential Vqq in the general metric background (4). To do this, as in Ref. [46], we should 
consider two kinds of string configurations: one is a connected configuration and the other 
is a disconnected one. 

For the connected configuration, the string starts at the AdS boundary r = 1/e and 
x = £/2, goes into the bulk, turns abound at r = and x = 0, and finally reaches at 
r = 1/e and x = —1/2. For the disconnected configuration, two straight strings start at 
r = 1/e with x = ±£/2, respectively and end at the tip of the background soliton, r = r . 
Here 1 /e is the UV cutoff and £ is the distance between two quarks. 

We now take the gauge t — r and x(r) = a. Then the induced metric Xab can be 
expressed as 

ds 2 = Xab dx a dx b = g^^-dx a dx b = i^r-A^) 2 + Ada 2 - r 2 f(r)dr 2 . (20) 

dx a dx" r 2 g(r) da 
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Thus, we need to minimize the following Nambu-Goto action functional 



S c x °- n t [r] 



drda* 



L 2 , dr 



r 2 g(r) da 



2 _|_ j»2jy.2yVj,j 



(21) 




Figure 8: The heavy quark potential as a function of distance I in the case ( = 0.2. In the 
right part of the plot, lines from the top to down are for /i < /i c (blue), fi ~ 2.134 (red), 
fi ~ 3.025 (black), fi ~ 4.920 (purple) and /i ~ 5.655 (green), respectively. 



The integrand can be regarded as a Lagrangian and a as time. Taking advantage of 
the fact that the Hamiltonian is conserved, we can deduce the equation of motion 



dr 
da 



" L V rtf(n) 



(22) 



where we have used the condition that the geometry is smooth at r(x = 0) = r*, i.e., 
dr / da\ r=r ^ = 0. Substituting (22) to (21), we obtain the on-shell world sheet action in the 
connected case 



ccon 
°x-t 



2T 
2ira' 



dr\ 



'fir) 



L 



9{r l_ / r 4 /(r) 



(23) 



where the UV cutoff e has been taken into consideration and w con is a finite part. The 
distance £ is related to r* by 



dr- 



L 



(24) 
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For the disconnected configuration, with using the gauge r — t and a = r and x'(r) = 0, 
we obtain the regularized on-shell world sheet action as follows. 



«~ = §5 1 drL M = + wd " n - <25) 

One can see that and S^ l 11° n have the same UV behavior. The divergence in (25) 
just manifests the infinite mass of heavy quarks. The divergence in (23) can be subtracted 
by considering the infinite mass of quarks in (25). And thus we can obtain the heavy 
quark potential w = w con — nj dtscon as a function of the distance I. We find that the 
behavior of the heavy quark potential is quite similar for different parameters, i.e., /3, ( 
and \i. We present Figure (8) with £ = 0.2 as a concrete example. One can see from the 
figure that in all cases, the heavy quark potential goes roughly linear with the distance 
in large I region, which implies the dual field theory is a confinement phase, while in 
the small I region, the potential shows a Coulomb potential behavior, which means that 
the dual field theory is in a deconfinement phase in the small i region. There is an 
associated confinement /deconfinement phase transition when one changes the distance £. 
This feature is consistent with the one shown in Figure (3) through the entanglement 
entropy calculation. Here the confinement /deconfinement transition point can be identified 
with the point where the potential vanishes. 



1.095 - 
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( 1 
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1 .080 { I : 

2 3 4 5 6 7 8 
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Figure 9: The critical distance l c as a function of chemical potential. Lines from right to left 
are for £ = (purple), ( = 0.2 (green), ( = 0.65 (red) and ( = 1.6 (blue), respectively. The 
dotted part of the lines for ( = 0.65 (red) and ( = 1.6 (blue) are not thermodynamically 
favored, trace the physical curves by always choosing the solid lines. 

In Figure (9) we plot the relation of the critical distance l c with respect to chemical 
potential. One can see that l c increases monotonically when ( is very small, then a peak 
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emerges and the position of the peak moves to the left as we increase £, finally, the peak 
is cut off for very large £ because there is a jump at fi c for the first order phase transition. 
One interesting point is that the critical distance is not always a monotonic function of 
chemical potential. We think that this arises due to the different behavior of the heavy 
quark potential in the small and large I regions. In Figure (10), we plot the heavy quark 
potentials versus chemical potential in the deconfinement phase and confinement phase, 
respectively. We can see clearly that the behavior of the heavy quark potential is quite 
different in two cases. 




23456 23456 



Figure 10: The heavy quark potential as a function of chemical potential at fixed £/L — 0.9 
(left plot) and £/L = 2 (right plot). Lines are for ( = (purple), ( = 0.2 (green), £ = 0.65 
(red) and ( = 1.6 (blue), respectively. The dotted part of the lines for ( = 0.65 (red) and 
C = 1.6 (blue) are not thermodynamically favored, trace the physical curves by always 
choosing the solid lines. 



5.2 Spatial Wilson Loop 

In terms of gauge/gravity duality, to calculate the VEV of a spatial Wilson loop, we 
consider a rectangular loop C along two spatial directions (x,y) [47] [48]. The expectation 
value of the loop can be obtained by the following AdS / CFT dictionary 

V, = J DXe~ s *-K (26) 

Here X denotes the series of world sheet fields and S x _ y stands for a world sheet action. 
Again, in the large N limit, the saddle point approximation is valid, and we only need 
to take the minimum value of the Euclidean action among the saddle points. Given the 
background (4), we can derive the expectation value of the spatial Wilson loop by using 
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the Nambu-Goto action for S x - y 



S x - y = J d 2 rj^DetXab 



dxdyJr 2 (r 2 + — , \ ), (27) 



2ira' J y r 2 g»(r) 

with a' the string tension and Xab the induced world sheet metric with a, fe the indices in 
the (rj 1 = x,rj 2 = y) coordinates on the world sheet. We take one of the spatial direction 
y goes to infinity. The quark and anti-quark are located at x = ±|, respectively. As the 
case with the temporal Wilson loop, in this case, there are also two configurations. For 
the connected configuration, we can make use of r'^f ~ L = —r 2 and write down equation 
of motion for r'(x) 



r 4 

„2 



r 2( r 2 _|_ L 2 r'{x) 2 



(2* 



Here we have used the following boundary condition r(x = 0) = r#, r'(x = 0) = 0. We put 
the solution of the equation of motion back to the world sheet action (27) and obtain 

s con = _2*1 f 1 j _L f_ = rEj^/i , co n) (2Q) 



where Y denotes the separation of quark and antiquark in y direction and the distance i 
in x direction is given by 

dr—^=^-A-z- (30) 

Similar to the case for the temporal Wilson loop, we can obtain a finite spatial heavy 
quark potential (pseudo potential) by subtracting the mass of two heavy quarks. The latter 
is related to the disconnected configuration. For the disconnected configuration, we take 
the gauge rj 1 = x and rj 2 = y and x'(r) = 0, and obtain the on-shell world sheet action 

SfT = 7^7 / £ dr-±= = + w d —)- (31) 

y 2na' J ro ^Jg{r) 2na' e 

We find that the behavior of the pseudo potential w = w con — yj dlscon with respect to I 
is very similar to the one in the case of temporal Wilson loop. We therefore do not show 
the behavior of the pseudo potential here. Instead, we plot in Figure (11) the critical 
distance l c with respect to chemical potential. Here the critical distance l c is defined by a 
vanishing pseudo potential as the case in the temporal Wilson loop. We can see that i c is 
non-monotonic and a peak appears even when ( is vanishing, which is a little bit different 
from the case in Figure (9). For large chemical potential, the critical distance increases 
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monotonically. Comparing with Figure (7), we find that the critical strip width £ cr u from 
entanglement entropy and the critical distance £ c from the spatial Wilson loop are much 
qualitatively similar with each other. This is expected since both quantities come from the 
spatial sector of the background soliton solutions. 




2 3 4 5 6 7 



Figure 11: The critical distance £ c as a function of chemical potential in the case of spatial 
Wilson loop. Lines from right to left are for £ = (purple), ( = 0.2 (green), ( = 0.65 
(red) and ( = 1.6 (blue), respectively. The dotted part of the lines for £ = 0.65 (red) and 
£ = 1.6 (blue) are not thermodynamically favored, trace the physical curves by always 
choosing the solid lines. Note that the curve for ( = (purple) will continue to increases 
for sufficiently large chemical potential as other curves. 

But we find that the behavior of the pseudo potential with the change of chemical 
potential is different from the case in the temporal Wilson loop. In Figure (12) we plot the 
pseudo potentials with respect to chemical potential for fixed distance £/L — 0.9 (in the 
deconfinement phase) and £/L = 2 (in the confinement phase). In the case oi£/L = 0.9, we 
see that as ( is not very large, the behavior of heavy quark potential is qualitatively similar 
to each other, no matter the order of the phase transition. At the beginning of the phase 
transition, it firstly decrease in the superconducting phase and then increases continuously 
for large chemical potential. The position of the bottom moves to small chemical potential 
side with the increase of (. For sufficient large (, the bottom is cut off and the pseudo 
potential increases monotonically with the increase of the chemical potential. 

6 Conclusion and Discussions 

To better understand the properties of holographic superconductor, we investigated the 
behaviors of two non-local physical quantities, i.e., entanglement entropy and Wilson loop 
in the Stiickelberg holographic insulator/superconductor model. The model exhibits rich 
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Figure 12: The pseudo potential as a function of chemical potential at fixed £/L — 0.9 
(left plot) and £/L = 2 (right plot). Lines are for ( = (purple), ( = 0.2 (green), ( = 0.65 
(red) and ( = 1.6 (blue), respectively. The dotted part of the lines for ( = 0.65 (red) and 
( = 1.6 (blue) are not thermodynamically favored, trace the physical curves by always 
choosing the solid lines. 

phase structure. For small (, the transition at fi c is second order and no discontinuity 
within the superconducting phase. For intermediate £, there is an additional complication. 
As we can see from the left plot of Figure (2), the grand potential as a function of chemical 
potential develops a "swallow tail" at a certain chemical potential /io in the superconducting 
phase. Thus, there is a jump in the condensate within the superconducting phase. For 
large £, the superconducting phase transition becomes first order. 

We calculated the entanglement entropy for a strip geometry and the VEV for temporal 
Wilson loop and spatial Wilson loop in the insulator /superconductor phase transition. 
Both show that there exists a confinement /deconfinement phase transition (see Figure (3) 
and Figure (8)). In both cases, we found that the critical width for the strip and the critical 
distance between quark and antiquark are non-monotonic functions of chemical potential 
(see Figure (7) and Figure (9)). 

There is either a discontinuity of the entanglement entropy or its slop with respect to 
chemical potential at /i c , namely, at the insulator/supercondcutor transition point, which 
indicates some kind of significant reorganization of the degrees of freedom as new degrees of 
freedom are expected to emerge in the superconducting phase. These discontinuities of the 
entanglement entropy or its slop correctly indicate the order of associated phase transition. 
Beyond the critical chemical potential, we found that the entanglement entropy is not 
monotonic in the superconducting phase: at the beginning of the transition, the entropy 
increases and reaches its maximum at a certain chemical potential [i m ax and then decreases 
monotonically. This behavior is universal in the holographic insulator/superconductor 
model. As one increases the model parameter (, but keeps the parameter /3, the position 
of the peak fi ma x moves to the small chemical potential side. 

Compared to the phenomenon observed in the entanglement entropy, the (pseudo) 
heavy quark potential with respect to chemical potential in superconducting phase presents 
more rich behaviors. For a fixed chemical potential, the heavy quark potential always 
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increases with the distance between the quark-antiquark pair (see Figure (8)). In the 
small i region the potential has a Coulomb potential form, while in the large £ region, 
the potential shows a confining potential form, as it is expected. But the behaviors of the 
heavy quark potential with respect to chemical potential are just opposite to each other in 
the small I region and in the large I region. 

In the spatial Wilson loop case, the pseudo potential shows a non-monotonic behavior 
with respect to chemical potential in the small i region, but not in the large I region (see 
Figure (12)). The non-monotonic behavior might be related to the one in the entanglement 
entropy. 

In this paper, we just reported the behaviors of the entanglement entropy and heavy 
quark potential in the holographic Stiickelberg insulator/superconductor model. Unfor- 
tunately, due to the lack of the knowledge of the dual field theory in the "bottom-up" 
approach, we can not completely understand some implications of the behaviors of entan- 
glement entropy and heavy quark potential, in particular, for the non-monotonic behavior, 
although we tried to give some discussions for this in Ref. [34]. But anyway, different from 
some physical quantities, such as condensate and charge density, in superconducting phase, 
we can see that there are various phase structures in superconducting phase by investigat- 
ing the behaviors of entanglement entropy and Wilson loop. It turns out that these two 
non-local physical quantities are good probes to study the properties of the holographic 
phase transition. 

In this paper we limited ourselves to the case with ^(i/j) = ip 2 + (tp 6 . For other 
form of the function, see, for example, Refs.[36, 38, 39], we expect that the result will 
be qualitatively same. But it is interesting to confirm this. In addition, It would be of 
great interest to investigate the behavior of entanglement entropy in the holographic p- 
wave insulator/superconductor case, in order to see whether the non-monotonic behavior 
is more universal or not. 
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